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Anderson localization (AL) is a ubiquitous interference phenomenon in which 
waves fail to propagate in a disordered medium. We observe three-dimensional 
AL of non-interacting ultracold matter by allowing a spin-polarized atomic 
Fermi gas to expand into a disordered potential. A two-component density 
distribution emerges consisting of an expanding mobile component and a non- 
diffusing localized component. We extract a mobility edge that increases with 
the disorder strength, whereas the thermally averaged localization length is 
shown to decrease with disorder strength and increase with particle energy. 
These measurements provide a benchmark for more sophisticated theories of 
AL. 

Wave propagation in disordered media is affected by scattering from random impurities. 
When those scattered waves self-interfere destructively, a phenomenon known as Anderson 
localization (AL) can arise (1). AL occurs in a wide variety of classical and quantum materials, 
impacting the transport of light (2, 3), acoustic (4), and matter (5, 6) waves. Localization is 
known to affect electrical conductivity in solids as a result of scattering from impurities and 
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defects (7), which is relevant to technological applications. We investigate 3D AL of a non- 
interacting ultracold atomic Fermi gas in a disordered potential created using optical speckle. 
The behavior of the gas is qualitatively consistent with several features of 3D AL and is shown 
to be incompatible with simple trapping and classical diffusion. Whereas in sufficiently large 
ID and 2D systems particles can be localized no matter how weak the disorder, AL is not 
inevitable in 3D (8). This distinctive feature of 3D AL — that the disorder strength sets a critical 
energy, the mobility edge, below which states are localized — is reflected here in measurements 
of the fraction of localized particles. We perform a direct measurement of how the mobility 
edge depends on disorder strength by smoothly tuning the speckle intensity. 

We create trapped, ultracold gases of fermionic 40 K atoms using standard techniques (9-11). 
The atoms are confined in an optical dipole trap, cooled to 170-1500 nK, and then spin polar- 
ized (12). Quantum statistics do not play a significant role in the measurements discussed here 
because the lowest temperature we sample corresponds to roughly one-half of the Fermi tem- 
perature of the gas. A focused optical speckle field (Fig. 1A) created by 532 nm light scattered 
through a diffuser and consisting of randomly distributed light and dark regions is superimposed 
on the trapped gas (as in (13)). The atoms experience a repulsive potential proportional to the 
speckle intensity, resulting in a disordered potential characterized by an approximately Gaus- 
sian autocorrelation function with ( x = 270 nm and ( z = 1600 nm r.m.s. radii. The laser beam 
creating the speckle propagates in the vertical z direction; we refer to the transverse directions 
(residing in the focal plane of the speckle) as x. The speckle intensity varies somewhat across 
the gas since it has a Gaussian envelope with a 170 /im 1/e 2 radius along x and y and a 400 /im 
Rayleigh range along z. The disorder strength A, which is the potential energy averaged over 
Cx and Cz at the center of the speckle field, can be continuously varied from 0-1000 k B xnK by 
adjusting the 532 nm laser power (k B is Boltzmann's constant). After the speckle field is slowly 
turned on over 200 ms, the trap is suddenly turned off, and the gas is allowed to expand in the 
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disordered potential while supported against gravity by a magnetic field gradient. Until the ex- 
pansion, the Gaussian momentum distribution of the trapped gas is unchanged by the presence 
of the speckle potential (12). 

Our experiments are distinguished primarily in two ways from previous work on AL of 
ultracold atoms. First, we work in 3D, where AL relies on small angle scattering rather than 
partial back reflections as in ID. Also, we eliminate inter-particle interactions by using a spin- 
polarized gas of fermionic atoms at temperatures far below the ~ 150 fiK p-wave collision 
threshold (14). In previous experiments, bosonic atoms were employed and the effects of in- 
teractions were suppressed or eliminated by using a Feshbach resonance (6) or by reducing the 
density (5). 

AL in 3D is conditional on the Ioffe-Regel criterion, which is equivalent to the quantum 
wavelength of the particle exceeding the Boltzmann transport mean-free path i B . For the maxi- 
mum A we achieve and the range of particle energies we sample, l B reaches a lower limit set by 
the speckle correlation length (15). Using £ B m (CxCzf^ and the thermal deBroglie wavelength 
A<iB — h/ \f2nmk B T, the Ioffe-Regel criterion corresponds to T < 300 nK. Here, h is Planck's 
constant, T is temperature, and m is the atomic mass. Since a spread of particle wavelengths 
are present in the gas and the Ioffe-Regel criterion is not a precise constraint, localization is 
possible even for temperatures somewhat above this limit. 

We probe localization by imaging the density profile after the gas has expanded for a vari- 
able time in the speckle potential. As observed for the typical absorption image (Fig. IB), a 
two-component profile emerges for any finite disorder energy. The mobile component has a 
Gaussian profile similar to that of a freely expanding gas, but expands (at a constant velocity) 
more rapidly than a thermal gas. In stark contrast, the stationary localized component has a 
profile along z that is approximately exponential, and a distribution along x that is well de- 
scribed by a Gaussian. Although strongly localized single particles are known to generally 
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possess exponential density profiles (7) with localization lengths £ that depend on the particle 
energy E, A, and the microscopic disorder parameters (15), a theoretical distribution applicable 
to our experiment (e.g., accounting for a thermal average over particle energies and localization 
lengths) is unresolved. Density profiles are therefore analyzed using a heuristic fit that repro- 
duces their basic features (12). The dynamics of the localized component are measured (after 
subtracting the mobile component from images (12)) using a fit to a distribution proportional 
to e^ 2 / 2 "^-! 2 !/^ with an exponential localization length £ z and r.m.s. radius a x along x (Fig. 
1C,D). 

As shown in Fig. 2, the size of the localized component becomes fixed after it rapidly 
expands along z for ~ 25 ms; the transverse size is approximately constant at the in-trap size. 
This apparent lack of diffusion cannot be explained classically. In a 3D speckle field there are no 
local intensity minima that can trap and classically localize particles, in contrast to ID (16, 17) 
and 2D (18, 19). Rather, a 3D speckle field consists of a rich collection of topological features, 
such as dark optical vortex rings and lines that do not propagate in a single direction but instead 
wander in all three directions on the length scales associated with ( x and ( z (20). While the 
vortex rings can trap particles in a finite volume, calculations of the percolation threshold for a 
3D speckle field establish that less than 0.2% of the particles are classically trapped for all of 
the data presented here (21). The expansion of the localized component is also inconsistent with 
classical dynamics. We numerically simulated classical trajectories in a 3D speckle potential 
for a thermal ensemble of particles under the conditions used for the data in Fig. 2C (12). The 
simulated sizes after expansion (solid lines) are incompatible with the observed dynamics of the 
localized component. 

Based on the expansion dynamics, we interpret the localized component as being comprised 
of particles with energies below the mobility edge E c , and atoms with higher energy constitut- 
ing the mobile component. The density profile of the localized component after expansion 
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can be straightforwardly understood in this context. The Gaussian distribution along x is sta- 
tionary because the transverse localization lengths are much smaller than the initial size of the 
gas; the observed profile cannot be explained by diffusion strongly suppressed via localization 
solely in z (12). The profile along z results from a thermal average of exponentially localized 
wavefunctions with much longer localization lengths, distributed such that the overall profile is 
approximately exponential. A disparity in localization lengths between x and z is expected be- 
cause i B , which controls the localization length and the range of energies over which it diverges 
when E ~ E c , strongly depends on the disordered potential correlation length. At low ener- 
gies, £b oc ( 2 , and at high energies £b oc ( 5 (in the weak scattering limit) (15), and therefore 
the typical localization length should be at least 36 times larger along z. 

By measuring the fraction of localized particles, we determine the mobility edge that sep- 
arates localized from extended states. The mobility edge we observe is unrelated to a similar 
quantity in ID (5) that results from correlations in the disordered potential (22). In ID, an effec- 
tive mobility edge arises because the quantum scattering generating back reflections becomes 
higher order above a momentum cutoff determined by the Fourier spectrum of the speckle po- 
tential. We measure the localized fraction after the atoms are held in the speckle potential for 
20 ms to resolve the mobile component and to minimize the impact of decay of atoms from the 
localized component (evident in Fig. 2A,B). This decay to finite localized fraction at long times 
is characterized by a 30-50 ms exponential time constant, does not strongly depend on tem- 
perature, and may in part arise from the atoms sampling lower intensity regions away from the 
center of the speckle field. Figure 3 A shows the fraction of localized atoms determined from fits 
to density profiles (12) for temperatures spanning 200-1500 nK and across the full range of ac- 
cessible disorder energy. More particles are localized as A is increased and E c correspondingly 
grows, or as T is decreased and fewer particles are thermally excited above E c . 

The mobility edge is determined from each point in Fig. 3 A by calculating the momen- 
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turn cutoff y/2mE c required to achieve the measured fraction of localized particles for a 3D 
spherically symmetric Gaussian momentum distribution consistent with the temperature (72). 
Figure 3B shows E c averaged across data taken at different temperatures for a fixed A. While 
E c increases with A, it does not follow the prediction that E c oc A 2 from the self-consistent 
Born approximation (23) and weak- scattering theory (15); for A > 80 /c^xnK the data fit well 
to a power law with E c oc A°- 59±0 02 (dashed line in Fig. 3B). Although these approaches are 
broadly applied to localization, this disagreement is likely a result of their failure precisely in 
the regime of AL. How A affects E c has not been predicted for the strongly localized regime 
that we probe here. 

The dependence of the measured localization length on the disorder strength and particle 
energy (controlled by changing temperature) is shown in Fig. 4. While data are only shown 
for a limited range of A and T, the observed qualitative behavior — that £ 2 increases with T 
and decreases with A — is characteristic of the entire range of parameters explored here. Ex- 
tracting universal parameters from the data shown in Fig. 4 is difficult in the absence of theory 
because the measured quantities are averaged across all particle energies present in the gas. The 
general monotonic trends, however, are consistent with a weak- scattering picture, in which the 
localization length is controlled by i B oc \fE / A 2 (15) at low energy. 

The variation of the speckle intensity across the gas is a complication that may affect the 
interpretation of the data shown in Figs. 3 and 4. While we have used the central speckle 
intensity A as a single parameter to characterize the disorder strength in Figs. 3 and 4, the 
speckle envelope causes a range of speckle intensities to be sampled by the atoms. The effect 
of the transverse speckle envelope is minor, since a x < 35 /im, which is small compared with 
the Gaussian waist. However, the maximum localization length in z is approximately 270 /im, 
which leads to the disorder energy averaged across the atomic density distribution reduced to 
approximately 70% of A; for smaller xi z the disorder energy averaged in this way is closer 
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to A. Although we have determined that this averaging does not directly affect the measured 
scaling of the mobility edge with disorder strength (12), the effect of the speckle envelope will 
likely be important to future comparisons between theory and our data. 

In the future, the exquisite control possible over ultra-cold disordered gases may enable 
measurements to shed new light on aspects of 3D localization that are not well understood or 
are complicated by inter-particle interactions or dissipation in other systems (see (24) for a re- 
view). An important issue that may be addressed if the single-particle states can be resolved 
(using Bragg spectroscopy, for example (25)) is the critical exponent that controls how the lo- 
calization length diverges for energies near the mobility edge. The influence of inter-particle 
interactions on the localization of fermions — another crucial question — can be studied by con- 
trollably introducing a second spin species; of particular interest is the impact of disorder on the 
BEC-BCS crossover (see (26, 27) and references therein). Finally, the impact of finite correla- 
tions in the disordered potential may be investigated using simple Gaussian or more complex 
holographic optics (28). 
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Figure 1: Ultracold gas expanding into an optical speckle field (green) and separating into 
localized (blue) and mobile (red) components. (B) The measured optical depth, proportional to 
the atomic density integrated through y, is shown in false color. The image is taken of a 480 nK 
gas that has expanded for 20 ms through the disordered potential with A = 240 fc^xnK. All 
images shown in this manuscript are averaged over at least 5 experimental realizations. Slices 
are shown through the image along x (C) and z (D). The filled curves are fits to independent 
mobile (red) and localized (blue) components. 
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Figure 2: Dynamics of the localized component for a 390 nK gas with A = 600 k B xnK. Slices 
through an image taken before release from the trap (black) and after the gas has expanded for 
40 ms (red) and 140 ms (blue) along the x (A) and z (B) directions. The wavelength of the 
imaging laser was changed to reduce the optical depth of the in-trap image by a factor of 15. 
The decrease in optical depth between 40 ms and 140 ms is a consequence of atoms slowly 
decaying from the localized component. (C) The measured localization length £ z (•) and r.m.s. 
size a x (■) of the localized component for variable hold time in the speckle potential. Each 
point is determined from an average of 6 images; the error bars (not visible for every point) in 
all figures are the standard error unless otherwise specified. The simulated classical expansion 
(solid lines) ignores rapid ballistic motion at short times, which lasts for several ms along z. 
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Figure 3: (A) The fraction of atoms in the localized component measured after 20 ms of 
expansion into the disordered potential for varying A and T = 240 ± 20 nK (•), 480 ± 20 nK 
(■), 1130 ± 60 nK (A), and 1470 ± 230 nK (♦). Each point is determined from fits to 5 images. 
The growing localized fraction with increasing A is evident in the insets to (B), which are 
images (with a false color logarithmic scale) taken at T = 480 nK and A = (i), 80 (ii), 160 
(iii), and 320 k B xnK (iv). (B) Using the data in (A), the mobility edge E c is determined at each 
A. Each point is a weighted average of E c accounting for the uncertainty in T and localized 
fraction. The error bars are the range of E c for the different temperatures that contribute to each 
point. 
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Figure 4: The dependence of the measured localization length £ z on disorder strength A at fixed 
temperature T = 480 nK ( ) and on T at fixed A = 480 /c^xnK (Q0- The green points are 
from the same data set as in Fig. 3 A, and the black points are an average of 10 experimental 
realizations. The error bars in T are from the uncertainty in the thermal expansion velocity. 
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Supporting Online Material 

Materials and Methods 

We prepare ultracold gases of 40 K atoms using techniques similar to those in (9) and (11). 
A mixture of atoms in the |F = 9/2, m F = 9/2) and |F = 9/2, m F = 7/2) hyperfine states 
is used for initial forced evaporative cooling in a variation of the QUIC magnetic trap (10). 
Transitions from both spin states to magnetically un-trapped Zeeman levels in the F = 7/2 
hyperfine state are driven using a single microwave-frequency magnetic field. After cooling in 
the QUIC trap, the atoms are transferred into a crossed-beam 1064 nm dipole trap. Following 
additional cooling induced by reducing the laser intensity, atoms in the \F — 9/2,m F — 7/2) 
state are removed using microwave transitions and a magnetic field gradient. The number of 
atoms in the gas ranges from 6 x 10 4 to 3 x 10 5 for all of the data in this manuscript, and 
the measured spin polarization is greater than 65:1. The waist of the Gaussian beam used 
for the dipole trap is approximately 120 /mi, and the trap frequencies vary from u z = u y = 
2nx 120 Hz, u x = 2n x 50 Hz to uj z = u y = 27rx210 Hz, u x — 2ir x 80 Hz during the final 
stage of cooling. The final temperature of the gas is correlated with the trap frequencies, with 
the lowest (highest) temperature corresponding to the lowest (highest) trap frequencies. 

The optical speckle field is generated by passing a 532 nm laser beam with a 16 mm waist 
through a 15 mm diameter plano-convex lens with a 13 mm focal length (Lightpath Industries 
GPX-15-15). A 0.25 mm thick holographic diffuser (from Luminit, LLC) mounted to the flat 
side of the focusing lens scatters the light through a 0.5 degree range of angles. The waist of 
the speckle field envelope is determined in situ by measuring the dipole force from the envelope 
on the mobile component, as in (13). We determine A within a 10% systematic uncertainty 
by calculating the average intensity based on the measured waist and 532 nm laser power. The 
Rayleigh length that characterizes the envelope along z is measured ex-situ. 



12 



Temperatures are determined by measuring the expansion velocity of the gas without dis- 
order present. The temperature and corresponding uncertainty are calculated from a linear fit 
to the r.m.s. size measured at different free expansion times. We determine the maximum im- 
pact of disorder on the temperature and momentum distribution using approximately the lowest 
temperature and highest disorder energy that we sample. We prepare a gas at (262 ± 4) nK in 
the dipole trap, then turn on the speckle potential, and then measure the momentum distribution 
after simultaneously turning off the speckle potential and the trap. For A ps 800 k B x nK, we 
determine that the temperature of the gas is (288 ± 9) nK. Thus, the speckle potential induces 
at most a 10% shift in the temperature. We also verify that the disorder does not affect the shape 
of the momentum distribution using a chi-squared analysis of the Gaussian fit to the time-of- 
flight images. We find that the reduced chi-squared is 1.26 ± 0.22 with A = k B x nK and 
1.32 ± 0.25 for A ps 800 k B x nK. The disorder therefore leads to no significant change in the 
shape of the momentum distribution. 

Images of the gas after expansion in the speckle potential are analyzed by first fitting only 
the mobile component to a Gaussian profile. This is accomplished by excluding the region 
containing the localized component using a rectangular mask. The size of the mask is approx- 
imately twice the size of the localized component along x and runs across the entire image in 
z. Small changes in the width of the mask do not affect the measured localized fraction, £ 2 , 
or a x . The mobile component is removed from the image by subtracting the Gaussian fit, and 
the residual localized component is fit to a function proportional to e^/^-N/^ The centers 
of the fits to the mobile and localized components are treated as independent free parameters. 
The localized fraction is determined from the fitted integrated column density in each compo- 
nent. Fringes that arise from technical noise are removed from the images shown in Figs. 1-3 
using standard image processing techniques; this procedure is not applied to images used for 
determining the properties of the localized component. 
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The simulation of classical trajectories is carried out for a range of particle energies, aver- 
aging over randomly sampled initial position. The computed dynamics are diffusive along all 
directions for expansion times greater than a few ms. To simulate the expansion of the gas, we 
average the diffusion constants over a thermal ensemble of energies consistent with the initial 
momentum distribution. The simulated r.m.s. radius along z is converted to a localization length 
for comparison to the data by assuming an exponential profile. The motion in x and y reaches 
the asymptotic regime within 10 /im in z. Therefore, for the localization lengths we measure, 
localization solely in z cannot explain the absence of diffusion in the transverse directions given 
the simulated diffusion rates. 

We determined that the speckle envelope does not directly affect the measured scaling of 
the mobility edge by analyzing the raw data from Fig. 3 using a different method. The disorder 
energy A averaged across the density profile is computed using the measured dimensions of the 
gas. For this alternative analysis, A is employed as a measure of the disorder strength instead 
of A. The mobility edge is computed for each point (at a specific temperature and speckle 
intensity) and plotted vs. A. The data are fitted to a power law, which gives E c oc A °- 58±0 03 5 
in agreement with the simpler analysis shown in Fig. 3B. 
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